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Abstract: We investigate the LHC sensitivity to supersymmetric models with light hig- 
gsinos, small R-parity breaking and gravitino dark matter. The limits on decaying grav- 
itino dark matter from gamma-ray searches with the Fermi-LAT put a lower bound on 
the higgsino-like neutralino NLSP decay length, giving rise to a displaced-vertex collider 
signature. Using publicly available tools for simulation of signal, background and detector 
response, we find that higgsinos with masses of 100 - 400 GeV and R-parity violation of 
( ~ 10~ 8 - 10~ 9 can show up in the 8TeV LHC data with 10 - 30fb _1 of integrated lumi- 
nosity. We demonstrate that in the case of a signal, the higgsino mass can be determined 
by reconstruction of the dimuon mass edge. 
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1 Introduction 

The Large Hadron Collider (LHC) at CERN has brought the field of high energy physics 
into a whole new era. The observation by both the ATLAS and CMS experiments of a 
new resonance at ~ 126 GeV [1, 2] is consistent with the long-sought Higgs boson. If this is 
indeed the scalar particle needed to induce electroweak symmetry breaking in the standard 
model (SM), the search for new physics does not end there. Such a light scalar receives 
quadratically divergent corrections to its mass, leading to large fine-tuning within the SM. 
Furthermore, although the gravitational evidence for dark matter is strong (see e.g. [3]), 
its particle nature has yet to be determined. 

The minimal supersymmetric standard model (MSSM) addresses both these issues but 
also faces new problems. While unbroken supersymmetry only introduces one new param- 
eter, the higgsino mass parameter fj,, soft supersymmetry breaking generically introduces 
many new parameters along with sources of large flavour- and CP-violation (see e.g. [4]). 



- 1 - 



The latter problems would be alleviated if the superpartners have large (multi-TeV) masses. 
Indeed, this would be consistent with a relatively large value of ~ 126 GeV for the Higgs 
mass and the absence of LHC signals for any other new particles so far. 

The disadvantage of such a decoupling solution is that heavy superpartners possibly 
reintroduce fine-tuning in the electroweak sector. Nevertheless, a light-higgsino scenario 
(LHS), in which two higgsino-like neutralinos and a higgsino-like chargino are light (of the 
order 100 GeV) and the other superparticles heavy, is a theoretical possibility. Experimental 
limits from chargino searches tell us that \i cannot be zero and there is a priori nothing 
relating its value to the soft breaking parameters. Even if the electroweak scale can be 
obtained with cancelations of terms with large soft parameters, this would be spoiled again 
if fi is too large (see e.g. [5]). 

In fact, a spectrum with higgsino masses around the electroweak scale and the other 
sparticles typically in the TeV range has been shown recently to arise in models of hybrid 
gauge-gravity mediated supersymmetry breaking [6]. Such models are motivated by ways 
of obtaining the MSSM and grand unification in string theory that feature a hidden sector 
of exotic states. They thus make use of the mechanism present in gravity- mediated su- 
persymmetry breaking to generate an electroweak-scale fx, as well as of the advantages of 
gauge mediation by providing a large number of messengers. 

In this set-up we have 

M « |Mi, 2 | , (1.1) 

where Mi 2 are the bino- and wino mass parameters, and mass splittings in the higgsino 
sector of the order m z/\M li2 \. Because the higgsinos are nearly mass degenerate and the 
strongly interacting superparticles are out of reach, such a scenario within the usual MSSM 
is difficult to probe at the LHC [7]. As will be the subject of this work, the prospects change 
if we allow for lepton number violating, and therefore R-parity violating, couplings. 

R-parity conserves baryon and lepton number and is imposed in the usual MSSM in 
order to forbid proton decay. Its conservation also renders the lightest supersymmetric 
particle stable, making the lightest neutralino a natural WIMP dark matter candidate. 
From a theoretical point of view, however, R-parity conservation is not particularly favoured 
and the stability of the proton can be ensured by requiring either baryon or lepton number 
violation to be small. In fact, it is a disadvantage of the MSSM LHS scenario that the relic 
density of neutralinos is too low, due to coannihilations, to provide the dark matter. Here, 
we will consider an extension of the MSSM with broken R-parity where the dark matter 
candidate is instead the gravitino. 

The possibility of gravitino dark matter has an interesting connection to leptogenesis. 
While the gravitino is a prediction of the desirable promotion of supersymmetry to a local 
symmetry, it leads to the gravitino problem [8-12]. The thermal production of gravitinos 
depends on the reheating temperature and the gravitino and gluino masses [13] . In thermal 
leptogenesis, the lepton asymmetry is created through the decays of heavy right-handed 
neutrinos and then transferred to a baryon asymmetry via sphaleron processes. In order 
to account for the small neutrino masses generated via the seesaw mechanism, the right- 
handed neutrinos need to be very heavy and therefore high reheating temperatures are 
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required to produce them [14-16]. This would also lead to gravitinos being produced in 
great abundance. If the gravitino is not the lightest supersymmetric particle (LSP), its 
decays would interfere with big bang nucleosynthesis (BBN) [17]. If the gravitino is the 
LSP and a dark matter candidate, the BBN bounds instead apply to the next-to-lightest 
supersymmetric particle (NLSP) which is meta-stable. This is where R-parity violation 
(RPV) comes in. 

The requirement of successful baryogenesis puts an upper bound on the amount of 
RPV that can be allowed, by the condition that the baryon asymmetry is not erased before 
the electroweak phase transition sets in. A small amount of RPV, however, leads to 1) 
a gravitino lifetime exceeding the age of the universe, because of the double suppression 
of gravitino decays by the Planck mass and the small RPV coupling and 2) a sufficiently 
short lifetime of the NLSP to be consistent with primordial nucleosynthesis. This makes 
it possible to have a good gravitino dark matter candidate even with the high reheating 
temperatures needed for leptogenesis, thus solving the gravitino problem [18]. 

It has been shown that the gravitino can account for the observed dark matter abun- 
dance for typical gluino masses and different types of NLSPs [19]. In the LHS, we will have 
to allow for a very large gluino mass and we discuss in this work how the gravitino dark 
matter abundance can be accounted for in this case. There are also direct bounds on the 
RPV couplings from cosmology, which will be of particular relevance for the present work. 
Decaying dark matter of this kind would lead to a diffuse gamma-ray flux observable by 
the Fermi-LAT. The non-observation of such an excess gives an upper bound on the RPV 
and thereby a lower bound on the NLSP decay length [20]. The finite NLSP decay length 
leads to the prediction of the displaced vertex signatures at the LHC that we study here. 

Displaced vertex signatures are extremely powerful in suppressing SM backgrounds and 
have been studied, in the context of different SUSY models, for macroscopic decay lengths 
ranging from 1 mm to hundreds of meters [21-24]. Unlike e.g. models where neutrino masses 
are generated by RPV, which lead to decay lengths up to 1 mm, the scenario described above 
predicts decay lengths that are orders of magnitude larger. This will lead to displacements 
of decay vertices in the outer layers of the multipurpose LHC detectors ATLAS and CMS, 
which motivates the muon signature that we consider. 

In this work, we consider the MSSM extended by bilinear R-parity violating couplings. 
These violate lepton number and R-parity by the introduction of only a small number of 
new parameters. Baryon number is conserved also at loop level, ensuring the stability of 
the proton. The bilinear interactions can be rotated under the symmetries of the theory 
to allow for a description in terms of trilinear Yukawa interactions, which simplifies the 
phenomenological analysis. In this framework, we study the case of a higgsino-like neu- 
tralino NLSP and the prospects for discovery with the data from proton collisions at a 
center-of-mass energy of 8 TeV. All the strongly interacting superpartners are assumed to 
be out of reach, and the higgsinos would be pair-produced via a virtual Z or W boson. 
The heavier neutralino and chargino will decay into the NLSP, and the NLSP will travel in 
the detector before decaying into SM particles, typically a W boson and a charged lepton 
in the case of a higgsino NLSP. NSLP decays into a gravitino LSP are, due to suppression 
by the Planck mass, orders of magnitude less probable. A clear signature arises from at 
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least one of the two NLSPs decaying inside the detector, giving rise to two opposite-sign 
muons from a secondary vertex. 

In section 2 we set up the framework of bilinear RPV, derive the branching ratios 
for the higgsino-like neutralino NLSP and the gravitino LSP and discuss the cosmological 
bounds on the RPV couplings. In section 3 we describe the LHC signature and analysis 
tools, and present the result of our detector level study for a few benchmark models. We 
conclude in section 4. 

2 Decaying dark matter in the light-higgsino scenario 

If we abandon the requirement of R-parity, the additional terms in the MSSM superpo- 
tential together with the soft terms introduce 99 new parameters into the model [25]. By 
allowing only for the bilinear terms, baryon number is conserved, and the number of new 
parameters is reduced to 9 (c./. [25, 26]). Such a scenario can be realised through the 
spontaneous breaking of B — L, the difference between baryon and lepton number [18]. 

Previous studies of RPV at the LHC have been motivated by the fact that bounds on 
sparticle masses can be weakened due to RPV, thus providing a possible explanation to 
why supersymmetry has not been discovered yet even if it is already being produced at the 
LHC [24, 27]. Here, we show that RPV can also allow for detection of weakly coupled new 
physics, which otherwise is not reached by the sensitivity of the LHC searches performed 
so far. Compared to the case studied in [24] where the lightest neutralino was assumed to 
be bino-like, direct production of higgsino-like neutralinos will have larger cross-sections 
since they are not suppressed by mixing angles [4] . In this section we derive the branching 
ratios of relevance for our phenomenological study and introduce the relevant parameters 
for the study of RPV in the LHS. 

2.1 Bilinear R-parity breaking 

In the extension of the MSSM with bilinear R-parity breaking that we consider, mass 
mixing terms between lepton and Higgs fields appear in the superpotential 1 , 

AW = tnHuLi , (2.1) 

as well as in the scalar potential, induced by supersymmetry breaking, 

-AC = BiHJi + m 2 J\H d + h.c. . (2.2) 

These mixing terms, together with the R-parity conserving superpotential, 

W = ^H u H d + hfjQiUjHu + hf 3 DiQ 3 H d + h^L^Hd , (2.3) 

the scalar mass terms, 

-Cm = m 2 u HlH u + ui\n\H d + (BH u H d + h.c.) 

+ ™>iiPji + mlicUi + mqiQiQi + fh 2 ui u\ui + m^djdj , (2.4) 

Our notation for Higgs and matter superfields, scalars and left-handed fermions reads: H u = (H u , h u ), 
Li = (li,h), Ei = (e, e), etc., where i is the family index. 
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and the standard SU(3) x SU(2) x U(1)y gauge interactions define the supersymmetric 
standard model with bilinear R-parity breaking. Note that the Higgs mass terms m\ and 
m\ contain the contributions both from the superpotential (2.3) and the soft supersym- 
metry breaking terms. For simplicity, we have assumed flavour diagonal mass matrices 
in (2.4). 

As discussed in [20], it is convenient to work in a basis of SU(2) doublets where the 
mass mixings /ij, B{ and mf d in Eqs. (2.1) and (2.2) are traded for R-parity breaking 
Yukawa couplings. This can be achieved by field redefinitions: the standard rotation of 
the superfields H d and li, 

H d = H' d - etl'i , k = l'i + uH' d , a = ^ , (2.5) 

A* 

followed by a non-supersymmetric rotation involving all scalar SU(2) doublets, 

H 'd = H d - e'il" , eH* = eH'* - e"l" , /■ = l" + t i H d + e'-eH'* , (2.6) 

where e is the usual SU(2) matrix, e = ia 2 , and e[ as well as e" are functions of B, Bi, m 2 d , 
ml, ml and m\ [20]. 

The R-parity breaking Yukawa terms contain couplings between gauginos, lepton dou- 
blets and Higgs doublets. After electroweak symmetry breaking one obtains new mass 
mixings between higgsinos, gauginos and leptons, 

-AC M ^ mtj—ejhZ - m z s w C*Vib + m z c w C*ViW Z + V2m z c w (*eiW + + h.c. , (2.7) 

where we have defined: 

efoi + efou V= J V 2 + V 2 ^ =tan/3 ^^ ; (2 .8a) 

v V Vd Cp 



rriij = hijVd , m z = -j= , s w = = = yjl - < . (2.8b) 

Here g, g' and are the SU(2) and the U(l)y gauge couplings and the charged lepton 
Yukawa couplings, respectively. Note that one also obtains couplings of the bino and wino 
to the lepton doublets and the Higgs doublets, 

-AC = -^= (e<flSV< + < H u^) b + ^= (4H° d *»i + <H°Ui) w 3 + h.c. , (2.9) 

where we have shown only the couplings to the neutral Higgs states. The neutral higgsinos, 
on the other hand, only couple to the charged Higgs. Introducing the physical Higgs fields 
in the unitary gauge, and taking only the coupling to the lightest Higgs into account, one 
obtains: 

—AC = —-gKihvib+-gKihviW + h.c. , (2-10) 

where 

Ki = e'jSin(— q) + e" cos(a). (2-11) 
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In the Higgs decoupling limit, a ~ (3 — 7r /2, we have Ki — Q. The Higgs decoupling limit is 
satisfied in the models considered in the present work. 

The details of the following derivation are given in Appendix A. The diagonal mass 
terms together with the mixing terms in Eq. (2.7) represent the 7x7 neutralino mass 
matrix of the gauginos b, u> 3 , the higgsinos h^, h° d and the three neutrinos v, L (A. 3a), and 
also the 5x5 chargino mass matrix of gaugino, higgsino and charged leptons (A. 3b). Both 
mass matrices have to be diagonalized by (bi-) unitary transformation matrices U in order 
to obtain the mass eigenstates of the neutralinos x° anci charginos > respectively. The 
currents which couple gauge fields to neutralinos and charginos are modified as well by 
these transformations and then depend on CKM-type matrix elements of neutral V^ x ° ,u \ 
charged V^ x °' e \ which are functions of the transformation matrices U. Furthermore we 
have derived the coupling to the supercurrents U^ x ° ,v \ as well as the coupling to the lightest 
Higgs V( u,x ). These R-parity breaking matrix elements are calculated in Appendix A and 
read 

= u% [ + w 2 ) ^ ~ <*> " ( Mi -,)(m 2 -^ + c ^ ) 

x [1 + J J , (2.12a) 



■/??, 



= [ 2 M 2 (M 2 -^ + <*> " (Mi — /x)(M 2 — |u) ^ + ^ " V, 



x 1 + s 2/3 ^f , (2.12b) 



Vi = AT 77—7 + ){sp + c p )\l + 0\ , (2.12c) 




^2 \M 2 -n M x -n, 
U^ ) =Qm z s w c w M ^(l + o{s 2 , r $\) . (2.12d) 




MiM 2 \ \ p m 2 

where the photino matrix element and the photino mass parameter are defined as 

C/p) = Cw U^ + s w U^ u) , = M lC ^ + M 2 4 , (2.13) 

and m is the largest out of the supersymmetric mass parameters Mi, M 2 and \x in the 
neutralino (chargino) mass matrix. Hence the approximated diagonalisation does not de- 
pend on the details of the supersymmetric spectrum, but in fact depends only on the ratio 
between the electroweak scale and the largest supersymmetric parameter. 

2.2 Gravitino and neutralino decays 

The partial width for a gravitino decaying into a photon and a neutrino (see Figure 1) is 
given by [28] 



2 rnl/ 2 
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Figure 1: Gravitino decay into photon and neutrino. 



Inserting the matrix element (2.12d) one obtains for the gravitino lifetime to leading order 
in m z/m [20] : 



r 3 / 2 (7i / ) 



(2.15) 



JL a 2^3/2 / M2-M1 V 

32 V / 2G F C M| ^ MiM 2 J ' 

where and a are the Fermi constant and the fine structure constant, respectively, and 
we have introduced the overall R-parity breaking parameter £ 



zZcf 

i 



(2.16) 



Contrary to the bino-like neutralino case (c./. [20, 24]) the gravitino lifetime cannot be 
expressed directly in terms of the mass of the lightest neutralino. Instead it depends on 
the mass scale of the gaugino mass parameter and hence on the masses of the heavier 
neutralinos. Nevertheless, one can invert the relation (2.15) with respect to the R-parity 
breaking parameter 



( = 4x 2 3 / 4 i 



IG F M P 



MiM 2 



a 



nr. 



M 2 - Mi 



(2.17) 



and constrain its value from the bounds on the gravitino lifetime using reasonable assump- 
tions for the bino and wino mass parameters (see Section 2.3.3). 

A neutralino with a mass larger than 100 GeV decays predominantly via two-body 
channels either into a W boson and a charged lepton, a Z boson and a neutrino [29] or a 
neutrino and a Higgs boson (see Figure 2). In the case of the higgsino-like neutralino the 
decay into a W boson and a charged lepton is dominant. For completeness we present the 
results for all three channels. The partial decay widths 2 read 

2 



r ( x ? -> w ± P 



Gf 3 



Zv 



hv 




1 + 2 



m 



w 



2 

Xi 



1 + 2 



m z 

2 

m 

Xi. 



a 



(2.18a) 



(2.18b) 



(2.18c) 



2 Decays into antiparticles have been taken into account. 
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Figure 2: Neutralino decays into neutrino and Z boson, charged lepton and W boson, 
and neutrino and the lightest Higgs boson. 

(y° e) (v° v) 

where V^* ' and ' are the charged and neutral current matrix elements, whereas 
Vh is the matrix element for the decay into the Higgs boson. We will evaluate these 
elements at the leading order using the expressions given in Eqs. (2.12a), (2.12b) and 
(2.12c), respectively. We will assume throughout the rest of this work that the Higgs 
decoupling limit is satisfied. 

Figure 3a shows the phase space suppression factors that are important for small 
neutralino masses. We have assumed that the mass of the lightest Higgs boson is 125 GeV. 
The total neutralino decay width is given by the sum 



r(xS -> w±F) + r( x ? -> zv) + r( x ? -> H 



(2.19) 



However, the evaluation of the partial decay widths in the case of a higgsino-like neutralino 
(see Figure 3b) leads to the conclusion that the lifetime of the lightest higgsino can be 
estimated solely from the decay into a W boson and a charged lepton. Below, we only 
provide the approximative formula, where we have taken into account that tan (3 in hybrid 
gauge-gravity mediation is in general large. In the numerical evaluation, however, we use 
the full results. Using Eq. (2.17) we can express £ in terms of the gravitino lifetime and 
arrive at the following expression for the higgsino-like neutralino lifetime 



.4 „4 



r x?- 



1 my 2 T 3/2C * s 



1 



(M 2 - Mi)' 



2tt 



Ml 



a 



fw(m x o) MlMlii 



m4 



M 2 (M 2 -n) 



0/3 + Cjs) 



My 



(Mi 



/i)(M 2 



(*j8 + c /?) 




(2.20) 



where /vk is the phase space suppression factor appearing in Eq. (2.18) 



fw(m x o, 



771 



W 



1 + 2 



m 




(2.21) 



and we have set m^o ~ ;U. The neutralino lifetime depends on the neutralino mass, the 
gravitino mass and its lifetime, and additionally on the bino and wino mass parameters. 
Expanding the higgsino lifetime in M/m allows us to arrive at a formula which is to first 
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(a) Phase space suppression factors for neu- 
tralino decays. 
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(b) Neutralino branching ratios. 



Figure 3: Phase space suppression factors and branching ratios for neutralino decays into 
W bosons (black), Z bosons (red) and the lightest Higgs (blue) and leptons. 



order independent of the higher neutralino mass scale m ~ M\ ~ M 2 : 



1 m- /2 n /2 si 6 

T — 



4 A 

w 



1 



1 



Ml 



M 2 - Mi 



M 2 sl 



l + O 



m Xi 



(2.22) 



a /H/(m x o)m x o \3Mi4 

The factor including the gaugino masses M\ and M 2 depends only on their ratio and is 
e.g. in the case of unified GUT masses very close to l /s. 

2.3 Cosmological bounds 

Both lower and upper bounds on £ can be derived from cosmology. The lower bound comes 
from the BBN contraints on the NLSP when the gravitino is the LSP. An upper bound can 
in principle be derived by demanding that the baryon asymmetry generated by leptogenesis 
is not washed out before the electroweak phase transition in the early universe. However, 
the bound from the constraint on decaying dark matter from the Fermi gamma-ray searches 
is stronger. As we will see, for our analysis the lower bound is not very constraining while 
the upper bound will be the motivation for our LHC search strategy. 

Having derived the decay widths of the gravitino LSP and the higgsino NLSP, we are 
now ready to estimate the gravitino mass range allowing for gravitino dark matter and 
successful leptogenesis [19]. This allows us to connect the results from gamma-ray searches 
with displaced neutralino decays at the LHC. 

2.3.1 Big bang nucleosynthesis 

To start with, we need to make sure that the decays of higgsino NLSPs do not interfere 
with BBN. Hence we demand that all higgsinos decay during the first 100 seconds of the 
universe [11, 30, 31]. Deriving the neutralino lifetime (2.22) without substituting the R- 
parity violating parameter £ for the gravitino mass and lifetime leads to 



C ~ 4.23 x 10 



-12 



,100 sj V 100 GeV 
which then characterizes the lower bound on £. 



X-3, 

2TeV 



(2.23) 
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2.3.2 Gravitino dark matter mass 

The minimal gravitino mass is limited by the requirement that the gravitino abundance 
does not overdose the universe. Since gravitinos are produced in thermal SQCD scatterings 
gg — > gfyfj, [32], the gravitino mass must increase with increasing gluino mass for a given 
reheating temperature. As we are interested in models in which the colored particles 
are inaccessible at the LHC, gluinos will typically be very heavy. For example, in the 
hybrid gauge-gravity mediation scenario in [33] which gives rise to a Higgs mass close to 
the tentative LHC result, the gluino mass is close to 4TeV. In order to still allow for 
small gravitino masses, we will assume that the hot phase of the universe was created in 
the decay of the false vacuum of unbroken B — L [34, 35] . Since right-handed neutrinos are 
created from B— L Higgs decays, this scenario allows for gravitino dark matter, leptogenesis 
and the correct values for the neutrino mass parameters while requiring lower reheating 
temperatures compared to the thermal leptogenesis case. The lower bound on the gravitino 
mass obtained in [35] for nig = 1 TeV is m™ n = 10 GeV. It is possible to scale this bound 
to other gluino masses using [34] 

»^=<(iw) a • (2 - 24) 

Assuming a lower bound on the gluino mass of 2 TeV, the minimal gravitino mass is 40 GeV 
and therefore a neutralino NLSP with a mass of 100 GeV is viable. 

2.3.3 Fermi-LAT bound on the gravitino lifetime 

With the help of Fermi-LAT data we are able to restrict the lifetime of gravitinos for a 
given mass. Using the isotropic diffuse gamma-ray flux one can derive a lower bound of 
T 3/2 ^ 3 x 10 28 s [20] . A stronger bound of r 3/2 > 6 x 10 28 s can be derived as a consequence of 
the non-observation of any gamma-ray lines [36-38] . For a LHS with a bino mass of roughly 
M\ ~ 2 TeV this translates via (2.17) into an upper bound on the R-parity violation of 
C < 4.70 x 10~ 8 and ( < 3.32 x 10" 8 , respectively. Finally we are able to derive a lower 
bound on the decay length of the lightest higgsino as a function of its mass as well as the 
mass and lifetime of the gravitino 

~* * — &r &) 3 (^?) ^k ? >-' . 

which is well within the reach of the multipurpose detectors at the LHC. Even when the 
higgsino decay length is larger than the detector dimensions, some higgsinos would, due to 
the statistical nature of the process, decay inside the detector. 

3 Higgsino production and decay at the LHC 

In the LHS, higgsino-like charginos and neutralinos would be pair produced at the LHC 
via virtual Z and W bosons. Heavier higgsinos decay into lighter ones, the lightest one 
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being the NLSP, typically a neutralino. The mass difference between the lightest chargino 
and the lightest neutralino is 

m z <m 1 2 / (1 + 828)(Mx — fj,) 2(? W (M 2 S 2 R + n) 

m+— m^o ~ —m 



* Xl - 2 z \(Mi -fi)(M 2 -fj) (M 2 + n)(M. 




M«m 1 2 / Mia + M 2 c 2 Mid + M|c 2 . r . , 
" M!M 2 + ^ M?Mf I . (3-1) 

where we have used 

ci = (1 - s 2/3 )<4 , c 2 = (1 + s 28 )sl . (3.2) 

The mass difference between the next-to-lightest neutralino and the lightest neutralino is 

m z <m 1 2 ( (1 + S 2 /3){My - fj) (1 - S 2 p) (My + jl) 



M<m 2 Msp / , /ll 1 



"Sis + ^ [m + m 2 ~ M,, ) + ° \k ) ) ■ < 3 - 3 ' 

Hence, in the case of heavy gauginos the mass difference is rather small and to first order 
proportional to m z/m. Therefore standard model products at this stage of the decay chain 
will be too soft to be detectable. In the presence of RPV, however, the NLSP would travel 
in the detector and further decay, prominently into a W boson and a lepton, yielding 
detectable SM objects coming from a displaced vertex. 

3.1 Signatures and search strategy 

As shown in section 2, neutralinos that are produced in proton collisions at the LHC decay 
in a secondary vertex into a W boson and a charged lepton in almost 100 % of all decays. 
Figure 4 shows an example of a cascade decay with muons in the final state. The distance 
between the collision point and the secondary vertex depends on the decay width of the 
neutralino (2.19), and hence on the R- parity breaking parameter £. 

The larger the value of £, the larger is the probability that the NLSP decays inside the 
detector. The different detector levels of the CMS detector are [39] : 

• The inner detector or tracker, for track reconstruction and vertex identification, which 
stretches out to a radius of r ~ 110 cm transverse to the beam. Its innermost part, 
the pixel detector, covers r < 11cm. 

• The electromagnetic calorimeter which measures electron and photon energies and 
stretches to r < 2 m. 

• The hadronic calorimeter, for measuring strongly interacting particles and identifica- 
tion of jets, which stretches out to r < 3 m. 

• The magnet, stretching out to r ~ 4 m. 

• The system of muon detectors or muon chambers, for identification of muons and 
measurement of their momentum, out to radius of r ~ 7.4 m. 
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invisible 




Figure 4: Typical R-parity violating decay chain involving higgsino-like neutralinos at 
the LHC. The secondary vertices as well as the two possibilities of interesting muon com- 
binations are highlighted. The Z boson decay is invisible, due to the small mass difference 
between the heavier higgsinos and the lightest higgsino (see Eqs. (3.1), (3.3) and Table 4). 
The signature is essentially the same for chargino production, since also in this case the 
decays into the lightest higgsino lead only to particles with small pt- 

Since the models under study here give rise to rather large displacements, we will choose 
to rely on the identification of muon objects in our search strategy. 

Table 1 summarizes all possible LHC signatures if the NLSP is a higgsino-like neu- 
tralino in the case when both neutralinos decay inside of the tracker volume. The signatures 
are classified according to the final states in the neutralino decays. Unlike the previously 
studied case of bilinear RPV in the MSSM in [24], when the spectra contain light squarks 
and gluinos and a bino-like neutralino NLSP, there are only two types of signatures: leptonic 
signatures involving only charged leptons in the final state, and semi-leptonic signatures 
involving at least two charged leptons and jets. 

Most supersymmetry searches for such final states at the LHC so far rely on the large 
production cross sections of the strongly interacting squarks and/or gluinos. The searches 
that are starting to probe direct electroweak production of charginos and neutralinos have 
been interpreted in R-parity conserving models with a stable neutralino LSP [40-42] . They 
rely on leptonic signatures with larger missing transverse energy than what would be 
present in our scenario in the case of both neutralinos decaying inside the tracker. Searches 
in the first LHC data for RPV as well as for long-lived particles have also been performed 
but because of differences in the scenarios considered or differing signatures they do not 
apply here. 

For smaller values of £, one of the neutralinos may decay inside or after the muon 
system leading to signatures with a larger amount of as shown in Table 2. We show in 
Figure 5 how this possibly gives rise to a missing energy signature as £ decreases. However, 
this situation would still not be covered by the current direct neutralino- and chargino 
searches by ATLAS and CMS [40-42] since the lepton reconstruction in these searches 
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category Xi decays LHC signature 



leptonic 


W+W~ 


i+r 


-> l+vl-vl+l- 








W+W^ 


-rr 


-> l+ul+vl'l- 


21- +21+ + j 






w~w~ 


1+1+ 


-> l~vl-vl+l+ 






semi-leptonic 


w+w~ 
w+w^ 


i+i- 
-rr 


-> 33l-vl + l- 
-> jjl+vl-l- 


2j + 2/- + 1/+ + J 






w+w~ 
w~w~ 


i+r 
-1+1+ 


-> jjz-w+r 


2j + 1/- + 2/+ + j 






w+w~ 


-ri+ 




4j + I/" + 1Z+ 




(same sign, no , 


$ T ) w+w " 

TJ w~w~ 


ri- 
-i+i+ 




4j + 2Z- 
4j + 2/+ 





Table 1: All possible final states in the higgsino-like neutralino case if both neutralinos 
decay inside the tracking volume. 

requires a track in the pixel detector, often with a certain maximum impact parameter to 
the primary vertex. For the same reason that we obtain missing energy from one of the 
neutralinos decaying outside of the detector, it will also be less probable that the other 
neutralino decays early enough for the leptons to fulfill such requirements on their inner 
tracks, as will be illustrated below. 

For very small R-parity violation both neutralinos may escape the detector without 
being observed. This case cannot be distinguished from the one with R-parity conservation 
and is for stable higgsinos and heavy coloured particles very challenging to detect [7] . This 
can be understood just by looking at the typical LHS process shown in Figure 4 and trying 
to imagine possible signatures in the case of the lightest neutralino being stable. 

We illustrate in Table 3 how often the different situations of none, one or two of the 
neutralinos, produced in processes like in Figure 4, decaying inside the detector occurs 
depending on the RPV and the higgsino mass. We also check how often the events would 
give inner tracks, possibly with a small impact parameter to the primary vertex, since this 
is a common requirement on lepton objects in existing supersymmetry searches. Table 3 
illustrates: 

• For decreasing C, we approach a situation indistinguishable from the case of R-parity 
conserving LHS (see numbers in red). 

• The second column in each subtable shows that when we can have large missing 
transverse energy in the event due to one of the neutralinos decaying outside of the 
detector, very few events will have the other neutralino decaying sufficiently close to 
the primary vertex for the event to be seen in existing searches for direct production 
of charginos and neutralinos (see numbers in blue) . 
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category 




X? decays 


LHC si§ 


nrature 


leptonic 
(opposite sign) 


W 
W 


n- - 
-i+ - 


-»■ l+ul- 


i/- + ir 


f + $ T 


single lepton 


w 


n- - 




2j + ir 


+ -^T 




w 


-i+ - 


-> ni + 


2j + ir 


f + $T 



Table 2 : All possible final states in the higgsino-like neutralino case if one of the neutralinos 
decays outside the tracking volume. 
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(a) For fj, = 100 GeV and ( = 1 x 10~ 7 , 
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(b) For u = 200 GeV and C = 5 x 10~ 8 , 
1 x 10~ 8 most often one neutralino leads to 
missing px- For smaller values of £ most 
of the time both neutralinos leave the de- 
tector undetected. 



Figure 5: Missing at the generator level, defined as the sum of the px of all neutrinos 
and of the neutralinos that decay after they have left the detector. In the case where only 
one of the two neutralinos escapes the detector, corresponding to a large value in the second 
columns in Table 3b, we notice an increased missing px around the neutralino mass. In 
the case where both neutralino decay outside the detector corresponding to a large value 
in the first columns in Table 3b, the values for small missing px are increased. This reflects 
that the neutralinos are mostly back-to-back, due to the absence of other high px objects 
in the decay cascade. 
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r < 5 mm 



Number of neutralino decays before the muon system 



1 2 

R-parity conserving LHS- displaced [J- + fJ>~ 
like 

1 may be covered by existing 

SUSY searches 

2 

Signature no LHC signature possibly large no or small $ T 

(a) Legend: Classification of event types depending on the number of neutralino decays inside 
the detector (columns) and inside the innermost 5mm of the tracker (rows). The fractions 
of events quoted in green belong to the type covered by our search strategy, the fraction of 
events in red cannot be distinguished from the R-parity conserving LHS. The fraction quoted 
in blue might be covered by existing SUSY searches (see text), and the black events might be 
misidentified as SM events. 
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(b) Generator level fractions in % of neutralinos decaying inside and outside of the detector 
as well as in its innermost part. The meaning of the positions in the subtables and the color 
code is explained in Legend (a). Numbers larger than 10% are bold. A zero entry means that 
< 0.001 % of the decays happen in this channel. The tendency of lighter higgsinos with smaller 
R-parity violation to decay outside the detector follows from relation (2.23). 



Table 3: Fractions of neutralino decays occuring either within a radius of 5 mm, inside the 
muon system or outside the detector, depending on the £ and n parameter. In Legend (a) 
we explain the color code and the meaning of the positions in Table (b) . 
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higgsino 








100 


200 300 


400 




106 


209 311 


413 


xf 


104 


207 309 


411 


X? 


102 


205 307 


408 



Table 4: Mass spectrum of light higgsinos in our four benchmark models with a higgsino 
mass parameter \x between 100 GeV and 400 GeV. All masses are given in units of GeV. 

• We also see that for a given £, increasing higgsino mass will lead to more decays 
inside the detector. 

The leptonic decays of at least one of the neutralinos inside the detector lead to an 
opposite-sign dimuon signature, which will be exploited in the present work. We will focus 
on events with two opposite sign muons originating either in a secondary vertex in the 
tracker, far away from the primary interaction point, or having no associated track at all, 
being identified solely by the muon system. The muon identification and reconstruction 
process applied is described in [24] and will be briefly summarized in Section 3.4. 

3.2 Benchmark points 

In models with hybrid gauge-gravity mediated supersymmetry breaking it is possible to re- 
alize a spectrum with higgsino masses around the electroweak scale, Higgs partners around 
the TeV scale and all other particles at the multi TeV scale [6]. 

The gravity mediated higgsino parameter of the superpotential and hence the light- 
est neutralinos and charginos can be of order 100 GeV. The existing lower bound on the 
chargino mass of roughly 95 GeV for degenerate spectra comes from LEP [43, 44]. In this 
study we have chosen four benchmark points for which we have varied the higgsino mass 
parameter \x in three steps from the lower bound of about 100 GeV to 400 GeV, see Table 4. 
The masses of the MSSM Higgs particles are in this case set by the CP-odd Higgs mass 
parameter m^, which we have taken to be 800 GeV. All other particles are governed by the 
gauge mediated parameters tuq and mu 2 which are chosen to be 3 TeV, putting them out 
of reach of the LHC. In all our benchmark points the Higgs mass is around 125 GeV. We 
have used these parameter choices as input values for a full RGE calculation performed 
with SDFTSUSY. As expected the production cross sections for all supersymmetric particles 
except the light higgsino states are negligible. The higgsino production cross sections for 
the four benchmark points are listed in Table 5. 

3.3 Background 

The SM processes that dominate the dimuon channel are: 
. tt 
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16.2 


a tot 


7649 
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25.73 


a tot 
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1371 
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18.3 


93.9 
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1940 




2.1 


2.6 


3.1 


3.7 


fgen 


565 


263 


1450 


5440 













(a) C = 8 TeV. (b) C = 14 TeV. 



Table 5: Partial and total NLO production cross sections for our benchmark models 
at 8 TeV and 14 TeV LHC in units of fb. The minimal and maximal (depending on Q) 
integrated luminosity corresponding to the generated number of events at 8 TeV for each 
model is given in units of fb -1 . 

. V*V*, where V = W,Z. 

In Table 6 we give the next-to-leading order (NLO) cross section for the processes that 
we have simulated for our study. As we will see, these will be efficiently removed by the 
requirement of a secondary vertex. 

In our analysis, we will require the muons to be isolated, which efficiently removes 
leptons originating in jets, and we further remove possible contributions from displaced b 
quarks by a sufficiently large cut on tracks in the inner radius. For low background levels, 
however, other sources of systematic uncertainties may come into play. These are: 

• cosmic muons, 

• pion and kaon decays in flight, 

• hadronic punch-throughs. 

An estimation of such contributions to our background has to be done with real LHC data, 
and is beyond the scope of this work. We argue here that most of this background, should 
it contribute, can be removed without significant loss of signal. Cosmic muons can be 
vetoed against by using the timing information, as discussed in [24], or a cut on back-to- 
back muons. Punch-throughs are also not simulated in Delphes but should in principle be 
possible to veto since in this case the muon would be associated with a jet. Most of any 
possible contribution to displaced muons from decay in flight should be removed by our 
high pt requirement on muons. In the following we therefore neglect the contribution of 
systematics to our displaced muon channel. We will, however, require our predicted signal 
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[45] 


[46] 


[47] 


[47] 


[47] 


CNLO 


183 


536 


57.25 


18.55 


7.92 


£gen 


196 


167 


360 


306 


143 



Table 6: NLO production cross sections for the relevant background processes in units of 
pb at an energy of 8 TeV as well as the integrated luminosity corresponding to the generated 
number of events in units of fb _1 . 

to amount to a certain number of observed events when our estimated background level is 
very low in order to allow some margin for systematics. 

3.4 Analysis 

3.4.1 Tools and settings 

All Monte Carlo samples were generated using MADGRAPH 4.4.44 [48] interfaced with PYTHIA 
6.4.22 [49]. We have chosen parton distribution functions given by CTEQ6L1 [50] and have 
used a modified version of SOFTSUSY 3.2.4 [51] for the calculation of the R-parity violating 
decays according to the formulas (2.18). All other signal decays were calculated using 
SDECAY [52]. 

The generic detector simulation DELPHES 1.9 [53], tuned to the CMS detector, was 
used in order to account for effects of event reconstruction at the detector level. The finite 
radial size of the detector important for studies with secondary vertices has been taken 
into account for the case of muons following the muon reconstruction process described 
in [24]. In the reconstruction process we have assumed that a muon can be reconstructed 
as long as its vertex is located before the muon system and that an inner track can be 
reconstructed as long as the vertex is in the first third of the tracker chamber. The gen- 
erator level particles and the detector level objects are matched including the distance in 
pseudorapidity and azimuthal angle. After the muon reconstruction process we are left 
with two non-overlapping classes of muon-objects: First what we will refer to as chamber 
muons, that are identified solely by the muon system, and second the tracker muons having 
in addition the information about their origin from the associated track 3 . The simulated 
muon reconstruction efficiency is close to 80 % in the inner parts of the detector and drops 
rapidly if the muon is created in the muon system. 

3.4.2 Cuts 

We focus our search solely on the muon objects as we assume that we can trust the detector 
simulation results in this case even in presence of secondary vertices. We aim to reconstruct 
clear signatures where one of the neutralinos decays into a muon and a W boson, which 
in turn decays into another muon and a neutrino. Thus, we demand two opposite sign 
muons that are either chamber muons or tracker muons with a secondary vertex far away 

3 In experimental jargon, these can also be referred to at standalone and combined muons, respectively. 
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Table 7: Outflow for the main SM dimuon background in units of fb. 



from the primary interaction point. We remind the reader that displaced muons are not 
expected in SM processes giving two isolated muons, and we will show in this section that 
the requirement of large displacements efficiently removes the SM dimuon background. 
Furthermore we demand that the invariant mass of the two-muon system is not too small, 
hence we suppress back-to-back signals. This reduces not only the SM background, but 
helps to decrease the background of cosmic muons. 

Only DELPHES muon objects with a px larger than 10 GeV are passed in our additional 
muon reconstruction processes taking into account the final radial size of the detector as 
described in the previous section. The cut on the pseudorapidity for muons is taken to be 
i] < 2.5. In addition to the isolation requirements in DELPHES, we only select muons that 
have no overlap with jet objects in angular coordinates AR > 0.1, where 

AR = ^A(f> 2 + At/ 2 , (3.4) 

and A(fr and Aij are the appropriate angle differences between the muon and the jet in the 
usual detector coordinates. Then we perform a selection cut on the total number of muons 
(chamber and tracker muons) in the event: 

• iV(muons) > 2. 

We define two event classes: 

• Class 1: the event contains exactly two tracker muons, i.e. iV(tracker muons) = 2, 
with opposite charge. 

• Class 2: the event does not fulfill the conditions for Class 1 and contains exactly two 
chamber muons, i.e. ^(chamber muons) = 2, of opposite charge. 

As the amount of R-parity breaking decreases, more and more events will fall into the second 
class. In accordance with the description of the signal presented above, we implement the 
following cuts on the Class 1 events: 
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1.64 


3.50 


0.82 


Class 1 


0.219 


0.109 


0.394 


0.096 


> 5GeV 


0.219 


0.109 


0.394 


0.096 


d(Vertex) > 5 mm 


0.219 


0.109 


0.394 


0.096 


Ad(Veitex)ij < 1 mm 


0.219 


0.109 


0.341 


0.075 


Class 2 


3.77 


1.31 


2.49 


0.692 


mirw{n + > 5GeV 


3.66 


1.31 


2.49 


0.692 


Total 


3.88 


1.42 


2.83 


0.767 



Table 8: Cutflow in units of fb for the two lighter benchmark points (/i = 100, 200 GeV) 
with the two most relevant of the analysed values of R-parity breaking £. 

• m irw{^ M~) > 5 GeV: We compute the invariant mass of the muon pair and demand 
that it is larger than 5 GeV. 

• d(Vertex) > 5 mm: Each of the tracks associated with the two tracker muons should 
have a vertex which is further than 5 mm away from the primary vertex. 

• A<i(Vertex)jj < 1 mm: The distance between the two track vertices should be less 
than 1 mm, in order to capture events where both muons originate in the same sec- 
ondary vertex. 

If the event fails one of the above cuts and fulfills the criteria for Class 2 events it is 
classified as Class 2: 

• "iinv(/U + M~) > 5 GeV: Also in this case we demand that the invariant mass of the 
muon pair is larger than 5 GeV. 

As expected, all of the SM background events are removed by the cut on the minimal 
distance of the vertex from the primary interaction point and the requirement that the 
reconstructed secondary vertices are close to each other (see Table 7). The LHS model 
events, however, survive these cuts to such an extent that a signal is detectable (see Table 8). 

3.4.3 Mass determination 

Rejecting the background-only hypothesis is only one contribution to the degree of belief 
that new physics has been discovered. Whether the signal hypothesis is a plausible one 
should be tested in other ways as well. Here, we show that the chosen signature allows for 
the determination of the neutralino mass via the well-known mass edge method [54, 55]. 
The mass edge in the dimuon invariant mass distribution is to leading order determined by 

ml = m 2 x o - m 2 w + O (j^J . (3.5) 
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Figure 6: Examples of the mass-edge reconstruction when the higgsino mass is m v o = 
102 GeV (benchmarks [i = 100 GeV), so that the theoretical value for the edge in the 
dimuon invariant mass is 62.4 GeV. The values in the legend show that already a small 
total number of events (#) are sufficient in order to reconstruct the mass edge with an 
error around 2 GeV. 



Following [56] we fold the phase space function with a Gaussian to model the mass edge: 
T(m u ) = -±= / o mCUt dy y exp (^) j , (3.6) 

where the endpoint m cu t and the height of the triangle a are the free parameters to be 
fitted to the dilepton invariant mass mu distribution to reconstruct the dimuon mass edge. 
We implemented this mass edge formula in the MINUIT class of the ROOT package. 

For this method to work, a sufficiently large sample of signal events is needed. In Fig- 
ure 6 we show examples of the mass edge reconstruction for different numbers of observed 
events in the case of our benchmarks model with \i = 100 GeV. We conclude from Figure 6b 
that a total number of events between 26 and 50 should give an accurate estimate of the 
higgsino mass. 

3.5 Result 

The signal cross sections after cuts for all our LHS benchmark models are given in Table 9. 
The SM background is removed by our cuts, as shown in Table 7. 

When dealing with very low background levels, a Gaussian approximation may not be 
adequate, and one should assume the number of events to be Poisson distributed. Under 
the null hypothesis of B background events, the probability of observing N or fewer events 
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c 



100 200 300 400 

5 x 10~ 8 90 ± 2.2 25.6 ± 0.52 4.94 ± 0.098 1.17 ± 0.025 

1 x 10" 8 3.88 ± 0.46 2.83 ± 0.17 1.07 ± 0.045 0.390 ± 0.014 

5 x 10~ 9 1.42 ± 0.27 0.767 ± 0.090 0.268 ± 0.023 0.105 ± 0.0073 

1 x 10~ 9 0.0283 ± 0.0071 0.023 ± 0.0093 0.0137 ± 0.0031 0.0037 ± 0.00082 



Table 9: Signal cross sections after cuts for all benchmark models and different values of 
the RPV parameter (, in units of fb. The errors are Poisson errors and the center-of-mass 
energy 8 TeV. 

is then 

P(N;B) = J2^ e~ B , (3.7) 

k=0 

given that the expectation value B is the true mean. We denote the expected number 
of events predicted by the model with S. To estimate the integrated luminosity 4 needed 
for a 5 sigma detection, one can require that there is a value for the minimum number 
of observed events iV bs = N + 1 such that 1 — P(N;B) < 2.9 x 10~ 5 %, corresponding 
to 5 standard deviations in the case of a one-sided Gaussian. In addition, -/V bs has to 
satisfy (1 — P(N; S + B)) being larger than some probability Pobs for observation under 
the hypothesis of our model. 5 

In our case, the expectation is B = 0, and in principle any S > would constitute 
a signal. In a real measurement, however, there will always be some background, and we 
will require N ^ s > 5 in order to have some margin to allow for systematic uncertainties. 
In Figure 7 we present results assuming a P ^, s of 50%, 90% and 99%. We see that 
the integrated luminosities of ~ 30 fb -1 expected in the 8 TeV run at the LHC would 
suffice for discovery of the lightest higgsinos with R-parity violation in the range C ~ 
2 x 10~ 9 - 2 x 10 -8 , and for £ above 6 x 10 -9 masses of \i = 400 GeV may be reached. 

After applying the mass edge method described in section 3.4.3 to one of our benchmark 
models, we estimate that 30 signal events are sufficient to reconstruct the neutralino mass 
with a couple of GeV's precision. Figure 8 shows the resulting integrated luminosities 
at which the higgsino mass could be determined at the LHC running at 8 TeV for our 
benchmark models. 

We use our 8 TeV results to estimate the reach when LHC runs at the design centre-of- 
mass energy of 14 TeV, by the same statistical analysis applied after a naive scaling of the 
cross sections after cuts with the factors "Ta/cj^ 1 presented in Table 5. Since the background 
is assumed to be completely removed by our cuts, the reach at 14 TeV would be significantly 
improved. Larger higgsino masses can be reached at smaller integrated luminosities, as can 

4 The integrated luminosity is C = aN where a is the cross section and N the number of events. 
5 For a discussion of the statistical measures used for this kind of study, see Appendix B of [57]. 
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Figure 7: Discovery reach with 8TeV center-of-mass energy at the LHC for our four 
benchmark models. Each colored band represents a value of [i and the lower, middle and 
upper line on each band corresponds to P Q b s = 50%, 90% and 99%, respectively. 




Integrated Luminosity L [fb ] 



Figure 8: Mass reconstruction reach at 8TeV under the assumption that S = 40 events 
(middle line of each colored band) are sufficient in order to reconstruct the neutralino 
mass. The lower and bands correspond to S = 30 and S = 50 events, respectively, and 
Pobs ~50%. 
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Figure 9: Estimation of the 14TeV discovery reach based on our 8TeV results. Each 
colored band represents a value of \i and the lower, middle and upper line on each band 
corresponds to P Q bs = 50%, 90% and 99%, respectively. 
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Figure 10: Estimation of the mass reconstruction reach for 14TeV based on our 8TeV 
analysis, under the assumption that S = 40 events are sufficient in order to reconstruct 
the neutralino mass. The bands correspond to the interval spanning between 30 and 50 
events, and P Q bs ~ 50 %. 
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be observed in Figure 9. Also the luminosity which is needed in order to reconstruct the 
neutralino masses is reduced, as can be seen in Figure 10. These results are approximate 
and we expect that changes in e.g. the px cut on the muons may be needed to deal with 
systematic effects such as increasing pile-up. 

4 Conclusion 

We have investigated the LHC detection prospects for the light- higgsino scenario, or LHS, in 
the MSSM extended with bilinear R-parity breaking terms. A spectrum with the higgsino- 
like neutralinos and chargino light, of the order of 100 GeV, and the other superparticles in 
the multi-TeV range is consistent with a Higgs mass of ~ 126 GeV within the MSSM and 
can be obtained in GUT models. Because the higgsinos are nearly mass degenerate and the 
strongly interacting superparticles are out of reach, such a scenario within the usual MSSM 
is difficult to probe at the LHC. The prospects change if we allow for R-parity violation, 
which leads to a consistent cosmology where leptogenesis and gravitino dark matter can 
be accounted for without conflict with BBN. 

We calculated the R-parity violating decay modes for a higgsino-like NLSP, and thereby 
diagonalized the full R-parity breaking neutralino matrix for this case. We showed that 
neither the resulting mass eigenstates nor the transformation matrices differ from the ones 
derived previously in the bino NSLP case. In fact, the parametrization only depends on a 
large hierarchy between the fermi scale and the largest neutralino mass. 

We were thus able to use the limits on decaying gravitino dark matter from gamma- 
ray searches with the Fermi-LAT to put an upper bound on the RPV parameter £ and 
thereby a lower bound on the higgsino NLSP decay length. This motivated our study of a 
dimuon LHC signature with macroscopic and, compared to other scenarios studied in the 
literature, relatively large displaced vertices. We simulated events and detector response 
for a few benchmark models, varying the value of the higgsino mass parameter \i and the 
RPV parameter £. We found that this scenario can show up already in the 8TeV run at 
the LHC with the 10 - 30 fb -1 of data expected at this center-of-mass energy, and that the 
reach in £ is improved by up to an order of magnitude compared to the current reach of 
gamma-ray searches. 

We also demonstrated that in the case of a signal, the LHS hypothesis can be tested 
further by a mass edge reconstruction. Except for the largest £ and smallest higgsino 
masses considered, the luminosity accumulated during the 8TeV run may not be enough. 
However, based on a our rough estimation of the 14 TeV reach, the higgsino mass in all our 
benchmark scenarios can be determined in the 14 TeV run, requiring integrated luminosities 
in the range 30 - 1000 fb" 1 . 

A Higgsino and gravitino branching ratios 

The couplings of the gauge fields to charged and neutral matter are given by 

C = -eJ^A" - ±J Zli Z» - -jftJ-W+f - -%J+W~» . (A.l) 



- 25 - 



The currents in the gauge eigenstate basis are 



T — T 3 4- T 2 ' 1 
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(A.2a) 
(A.2b) 

(A.2c) 
(A.2d) 



The upper indices of the electromagnetic currents indicate the transformation properties 
of the fields in the current under the SU{2)\ J . Having derived the currents we now have 
to transform them into the mass-eigenstate basis of the fermions. To this end, we have to 
diagonalize the mass matrices A4 N and A4 C , which in the gauge eigenbasis are 6 
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(A.3b) 



The matrices Ai N and A4 C are diagonalized by unitary and bi-unitary transformations, 
respectively, 



(A.4) 



where [/wt {/("■) = U^U^ = {/( c )t[/( c ) = 1. These unitary transformations relate the 
neutral and charged gauge eigenstates to the mass eigenstates (%°, z/{) (o = 1, ... ,4) and 
(Xa,e-), (xi,e'f) (a = 1,2), respectively, 
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5 Note the extra factors of y/2 in the charged mass matrix compared to [20] . 
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Our convention for the separation of the transformation matrices into R-parity conserving 
and violating parts is 
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(A.6a) 
(A.6b) 



After the transformation the currents can be expressed through CKM-like products of two 
transformation matrices 



r(x°)„o 



(A. 7a) 
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where the CKM-like matrices relevant for this study are given by 
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The following results are approximated by an expansion in £ and e, where e = m z/m and m 
the largest mass parameter of either Mi, M 2 or [i. The parameter choice affects neither the 
expansion of [/W nor the mass eigenstates. For the R-parity violating part of the neutral 
CKM-like matrix we find 
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Numerically, the relative errors are smaller than 0.10, 0.20, 0.15, 0.05 for a 
the R-parity violating part of the charged CKM-like matrix we find 
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Here we again neglected corrections that involve the Yukawa couplings h^. The numerical 
corrections to the NLO contributions to V^f'^ are smaller than 0.05, 0.15, 0.20 for a = 1, 
2,3, respectively. For a = 4 we reach the limit of our numerical precision. 
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